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Spontaneous breaking of chiral symmetry is investigated in the Hamiltonian approach to QCD 
in Coulomb gauge. The quark wave functional is determined by the variational principle using 
an ansatz which goes beyond the commonly used BCS-type of wave functionals and includes the 
coupling of the quarks to the transversal spatial gluons. Using the lattice gluon propagator as input 
it is shown that the low energy chiral properties of the quarks, like the quark condensate and the 
constituent quark mass, are substantially increased by the coupling of the quarks to the spatial 
gluons. Our results compare favourably with the phenomenological values. 
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Introduction. — The infrared sector of QCD is charac- 
terized by two non-perturbative phenomena: confine- 
ment and spontaneous breaking of chiral symmetry. For 
Nf massless quark flavours the QCD Lagrangian is in- 
variant under separate global flavour rotations of the left- 
and right-handed quarks 

UL{Nf) X Ur{Nj) = Uv{Nj) X Ua{Nj) . (1) 

The vector {V) and axial-vector (A) groups rotate left- 
and right-handed quarks in the, respectively, same and 
opposite way. In the quantum theory, UA^Nf) is sponta- 
neously broken by quark condensation {qq) ^ resulting 
in the generation of a constituent quark mass and the 
occurence of N'i Goldstone bosons. The latter can be 
identified with the (light) pseudoscalar mesons. Chiral 
symmetry is a good starting point for the Nf = 3 light 
quark fiavours u,d,s. The finite (current) quark masses 
explicitly break the axial group UA{Nf) and induce finite 
masses for the pseudoscalar mesons. 

Originally the dynamical breaking of chiral symmetry 
was studied within QCD-inspired models which share 
with QCD the property of chiral symmetry, e.g., the 
Nambu-Jona-Lasinio model |l|-|^ (which, however, is not 
confiniiig), or models with a confining two-body interac- 
tion [5h7] . In recent years there have been also attempts 
to study chiral symmetry breaking in continuum QCD 
either by means of Dyson- Schwinger equations [8|, |9| or 
by the renormalization group flow equations 'lO|. In this 
letter we study the spontaneous breaking of chiral sym- 
metry within a variational approach, which is the exten- 
sion to full QCD of the Hamiltonian approach to Yang- 
Mills theory in Coulomb gauge developed in recent years 

Variational calculations in Coulomb gauge Yang-Mills 
theory were started in Ref. llSlI and later on taken up 
in Ref. |16| . Our approach |lll4l4 1 differs from that of 
Refs. |15l . [iGJ in the ansatz for the vacuum wave func- 
tionals and, more importantly, in the full inclusion of the 
Faddeev-Popov determinant and in the renormalization, 
for more details see Ref. [l7|. Our approach has given 
a quite satisfactory description of the essential infrared 
properties of Yang-Mills theory: a linearly rising static 
quark potential [l3|, a gluon propagator in accord with 



Gribov's formula and with lattice data 18[, a perimeter 
law in the 't Hooft loop I4y and a dielectric function of 
the Yang-Mills vacuum [19| embodying the phenomeno- 
logical picture of the bag model. Similar infrared prop- 
erties were also obtained recently in a renormalization 
group flow equation approach to Hamiltonian Yang-Mills 
theory in Coulomb gauge [20|. 

In the present letter we extend the variational ap- 
proach in Coulomb gauge to full QCD. We use an ansatz 
for the vacuum wave functional of the quark sector which 
goes beyond the BCS-type wave functional used e.g. in 
Refs. [5|,|7| and includes the coupling of the quarks to the 
transversal spatial gluons [2l|. We will show that, using 
either the lattice results or the results of our variational 
approach for the gluon sector, this coupling changes the 
chiral properties of the quarks, like the quark conden- 
sate and constitutent mass, by 20-60% and brings these 
quantities in the region of their phenomenological values. 

Hamiltonian approach to QCD in Coulomb gauge. — 
The part of the QCD Hamiltonian in Coulomb gauge 
which contains the quark fields 'ip{x), is given by 24| 



H = H^ + H, 



QGC 



Ha 



with 



H„ 



i / d'^xilj'{x)[a.-d + l3mo]'il;{x) 



Hqcc = — g d'^xip''{x)a ■ A{x)^{x) , 
i/c = y f d^xd^yp^x)F^\x,y)p\y) 



(2) 



(3) 
(4) 



^ / d^xd^yJ-'p'^^(x)JF^'>{x,y)p\y) 



+ '- J d'xd'yp^{x)F'^\x, y)p^y„(y) . (5) 

Here a and /3 are the usual Dirac matrices, mg denotes 
the current quark mass, and A{x) = A°(a;)T° is the 
transversal spatial gauge field in the fundamental repre- 
sentation T" of the gauge group. Furthermore, 



p"(a;) =7/;t(a;)rXa;) 
p3v„(^) = r'''AUx)IlUx) 



(6) 

(7) 



are the quark and gluon color charge densities, respec- 
tively, with /"^'^ being the structure constants of the 
gauge group and Il'j.{x) ~ S/i6A'^{x) being the canon- 
ical momentum operator of the gauge field. Finally, 
J{A) — Det{—Dd) is the Fadeev-Popov determinant and 

F-\x,y) = {x\[{-Dd)-^ {-d^) {^Dd)-^Y'\y) (8) 

is the so-called Coulomb kernel where D'jf = S'^^dk + 
gj:acb^c jg ^-j^g covariaut derivative in the adjoint repre- 
sentation. 

The Coulomb term He, Eq. ([5|), arises from the resolu- 
tion of Gauss's law, so that gauge invariance is fully taken 
into account. The Yang-Mills vacuum expectation value 
of the Coulomb kernel 



{F-'{x,y)),.,,^5-'F{x,y 



(9) 



defines the static non-Abelian color Coulomb potential 
F{x,y). Lattice measurements show that this quantity 
rises linearly at large distances, however, with a coeffi- 
cient Cc (referred to as Coulomb string tension) which 
is about two to three times the value of the Wilsonian 
string tension Uw 25l-l28l|. In momentum space 



^F(fc) 



Sttctc 







(10) 



Due to asymptotic freedom at small distances the static 
non-Abelian color Coulomb potential F{x,y) behaves 
essentially like an ordinary Coulomb potential (up to 
anomalous dimensions). 

Variational approach to the quark sector of QCD. — 
We assume the QCD vacuum wave functional to be of 
the form |<i>) = \<j)) (g) |YM) with a Gaussian type-of wave 
functional for the Yang-Mills sector ll|, \l^ 



{A\YM) = J-Ha) exp (^-^ J At{xMx, y)A^{y) 



(11) 



For the quark sector we use the following trial ansatz 



: A/'exp 



ilMx)K{x,y)il}^{y) 



(12) 



where 'il}±{x) are the positive (negative) energy compo- 
nents of the quark field ip{x) and |0) is the vacuum wave 
functional of the free quarks defined by 



V'+(a;)|0) = = 7/'L(a;)|0) 



(13) 



With arbitrary kernel K the state | (f) is the most general 
Slater determinant which is not orthogonal to the vacuum 
|0). For the kernel K we assume the form 



K{x, y) = Ko{x, y)+ I dPz K{x, y; z) ■ A{z) , 



(14) 



where Kq and /C are determined by minimizing the en- 
ergy density of the quarks. Finally, A/" is a normalization 
constant, which is a functional of the kernel K, Eq. ()14|) . 
and thus of the gauge field. Through the vector ker- 
nel /C, the wave functional ([T^ contains the coupling of 
the quarks to the gluons. If this coupling is neglected 
(by putting K = 0) the wave functional \(j)), Eq. (IT4l) . is 
of the BCS-type used in Refs. (a LZ|, for which the ex- 
pectation value of the quark-gluon coupling term Hqqc, 
Eq. (jlj, vanishes. 

Our strategy is as follows: We first take the expecta- 
tion value of the quark Hamiltonian H, Eq. ([2]), in the 
fermionic state |0), Eq. p2|) . and subsequently in the 
Yang-Mills vacuum state (fTTj) . The fermionic expecta- 
tion value {H)p = {(j)\H\(p) can be taken exactly (by using 
Wick's theorem) since the wave functional \(j>) is a Slater 
determinant. It can be expressed in terms of the static 
quark propagator (in the presence of the gauge field) 



Gix,x')^^{^\[^Pix),4'Hx')] 



(15) 



which for the state (J12p is given in a matrix notation by 

G = A+(l + KK^)-^A+ + A_(l + K''K)-^K^A+ + 

+ A+(l + KK'^y^KA^ - A_(l + K^K)-'^A^ + 

A+), (16) 



1(A- 



where A± are the projectors on positive and negative en- 
ergy states V'± = A±i/j. Note that the quark propagator 
(and accordingly (i?)F) is via the kernel K, Eq. (fT4|) . a 
non-local and non-linear functional of the spatial gauge 
field A. 

For the gluonic expectation value (. . .)ym = 
(YM| . . . |YM) of {H}p in the Yang-Mills vacuum state 
([TT]) . we will restrict ourselves to two (overlapping) loops. 
To this order, in the gluonic expectation value of {H)p 
we can use the replacement 

(. . . (1 + K^K)-^ . . . )vM ^ (. . . (1 + (i^^if )ym)~' . . . )yM 

(17) 
and furthermore, replace the Coulomb kernel ^ by its 
vacuum expectation value, i.e. by the static non-Abelian 
color Coulomb potential ^. It is then straightforward, 
although quite involved, to evaluate {{H)f)ym and take 
the variation with respect to the kernels Kq^K. and lu. 
In the present study we ignore the back-reaction of the 
gluon sector to the presence of the quarks, i.e. we keep the 
gluonic kernel (the gluon quasi-particle energy) w at its 
form determined in the pure Yang-Mills sector |lll4l3j . 
We also ignore the small current quark mass rriQ, which 
can be easily included if desired. Furthermore, inspired 
by the form of the quark Hamiltonian H, Eq. ([2]), we 
assume for the quark kernels p^ of the wave functional 
(dH) the form 

Koix,y)=pSix~y) 
K.{x,y;z) = aV{x-y,z-y) , (18) 
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FIG. 1. The variational kernels S and V for ac = 2(Tw 

where S{x—y) and V{x—y, z—y) are scalar functions. In 
principle, more complex tensor structures of the quark- 
gluon coupling can be included as has been done in three- 
dimensional QED in Landau gauge j29| . 

Variation of the energy ((H)f)ym with respect to the 
kernels Eq. P^ 5, V yields the following gap equations 
in momentum space 



S{k) 

V{k,p) 

where 
R{k) = Cf 



i^Hk) 



|fe|-Vc(fc) 

g 1 + S'^{k) + R{k) 
2 |fc|-V,c(fe)+/F(fc) ' 



(19) 
(20) 



dh 



{2ny 



■D{£)V^{k,q) \ + {k-l){q 



with I = kg, Cf = {N^-l)/{2Nc) being the quadratic 
Casimir and D{£) being the static gluon propagator. In 
the above equations we have introduced the loop integrals 



Iqccik) ^ gC'F 



^ Di£)Vik,q)\l + ik ■£){<! ■ 



(22) 



9 



Cf 



d\ 



F{k 



2" J (2^fl + S^{q)+R{q) 
S{q){l-S^{k) + R{k))+ 

-{k-q)S{k){\-S\q)-R{q)) 



(23) 



li'\k) 



-Cf 



d\ 



F{k 



{2^fl + S^{q)+R{q) 



2S{k)S{q) + -{k • q) (1 - S^q) - R{q)) 



(24) 



FIG. 2. The dynamical mass M for the full solution and for 
V^ = with ac = 2(Tw. 



A closer inspection of Eq. ([^0]) shows that the vector ker- 
nel V{k,p) = V{k) does not depend on the momentum 
p which is conjugate to the coordinate z — y in (|18p . 
After taking the gluonic expectation value of the quark 
propagator (I16p one finds the quark condensate 



(^^) 



2Nr 



dqq 



S{q) 



1 + S^{q)+Riq) 



and the dynamical quark mass 



M(fc) 



2S{k) 



1 - S'2(fc) - i?(fc) 



(25) 



(26) 



Numerical Results. — For the static gluon propagator 
D{q) = l/{2uj{q)) we can take either the result of the 



(21) variational calculation [ll| - |l3l | or the lattice result [1 



Both propagators are very similar to each other and 
can be well fitted by Gribov's formula [3(| cj(fc) = 
yFTMVF with M « 880 MeV. 

The quark Hamiltonian H, Eq. ^, like the full QCD 
Hamiltonian, is scale invariant. In the present case the 
scale is set by the Coulomb string tension CTc, which en- 
ters the static quark potential F, Eq. PU)) . Further- 
more, the Hamiltonian H contains the coupling con- 
stant g explicitly. In principle, after full renormaliza- 
tion the coupling constant should enter only in the form 
of the running one. Since we are interested here in 
the infrared sector of QCD we will replace the run- 
ning coupling calculated in the Hamiltonian approach 
in Ref. 3l| by its (finite) infrared limit. Furthermore, 
we use aw = (440MeV)2 and iVc = 3. With this input 
the coupled equations (fT9|) and ([20|) are solved numeri- 
cally and the results are shown in Figs. [1] and [2] The 
quark condensate (|25l) and the dynamical mass (pS)) are 
affected in two ways by the inclusion of the coupling of 



the quarks to the spatial gluons: first through the change 
of S{k) and secondly through the occurence of the vec- 
tor form factor V{k) in the coupled equations P^ and 
(PU)) . Neglecting the coupling of the quarks to the gluons 
(/C = 0) we obtain for the quark condensate 



W) 



-nSMeVy/ac/a^ 



(27) 



which agrees with the findings of Ref. 7] if the same value 
for the string tension is used. With the inclusion of the 
quark-gluon coupling this value is shifted to 

(iPiP) « (-135MeVv/ac/aw)^ , (28) 

which is a 20% increase of the figure in the bracket. Lat- 
tice calculations show that CTc = (2 . . . 3)crw, which yields 
for the quark condensate ^E\\ 

(^V)«(-(191...234)MeV)\ (29) 

which compares favourably with the phenomenological 
value {'ipip) — (—230 Me V) . An even larger effect of 
the quark-gluon coupling is obtained for the constituent 
mass. Neglecting the quark-gluon coupling (/C = 0) 
yields 



M«84MeV\/CTc/o-w, 
while with this coupling included one finds 
M 



(30) 



132MeVv/crc/crw, (31) 

which is an increase of 57%. Using again the values for 
CTc quoted above we obtain 



Af w (186...230)MeV. 



(32) 



Conclusions. — Our results show that the coupling of the 
quarks to the spatial gluons, neglected in previous inves- 
tigations [5ri7|] , quite substantially influences the infrared 
chiral properties of the quarks. Typically the effect of the 
quark-gluon coupling is in the range between 20-60%. 
To make more accurate predictions within our approach, 
a more precise lattice measurement of the Coulomb string 
tension ctc is required, which sets the scale in our ap- 
proach. 
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